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Abstract 

We study the spectrum of the Dirac Hamiltonian in the presence of high temperature 
sphaleron-like fluctuations of the electroweak gauge and Higgs fields, relevant for the 
conditions prevailing in the early universe. The fluctuations are created by numerical 
lattice simulations. It is shown that a change in Chern-Simons number by one unit 
is accompanied by eigenvalues crossing zero and a change of sign of the generalized 
chirality T 5 = (— 1) 2T+1 7 5 which labels these modes. This provides further evidence 
that the sphaleron-like configurations observed in lattice simulations may be viewed 
as representing continuum configurations. 
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1 Introduction 

The baryonic current is anomalous in the electroweak theory 0. This leads to a 
non- conservation of baryon- and lepton-number. It is of little importance in the 
present day universe where the processes which violate baryon- and lepton-number 
conservation are exponentially suppressed. There is evidence that this suppression 
is not present at the high temperatures prevailing in the early universe f2j, and it 
has played an increasing role in attempts to explain the baryon asymmetry observed 
in the universe today. 

Unfortunately almost any interesting question related to the high temperature 
phase of the electroweak theory is difficult to address by analytical tools, due to 
the infrared singularities of the associated high temperature perturbation theory. 
A generic non-perturbative way to study the high temperature fluctuations of the 
gauge and Higgs fields is by means of lattice simulations. If we want to address ques- 
tions related to the anomaly we encounter the difficulty that topology is involved, 
an inherent continuum concept, and it is important to verify that configurations 
generated on the lattice during the simulations will qualify as representatives of 
continuum configurations. 

In [fj] a real time evolution of classical gauge and Higgs fields was used as an ap- 
proximation to the high temperature fluctuations present in the electroweak theory 
in the early universe. The baryon number violation was studied by observing the 
change in the Chern-Simons number during the time evolution. Since the temper- 
ature fluctuations are not small in these simulations it is important to ask if these 
fluctuations might spoil the picture of level crossing, which somehow is the justifi- 
cation of viewing the change in Chern-Simons number as responsible for a change 
in fermion number. 

In || it was shown that the gauge field fluctuations which lead to a change in 
Chern-Simons numbers at the same time cause the lowest eigenvalue of the massless 
Dirac Hamiltonian to cross zero, in agreement with the level crossing picture where 
such a crossing results in the conversion of particles to anti-particles (or vice versa). 
In addition the measurement of the spatial extension of energy lumps associated 
with the fluctuations corroborated the interpretation of these as representing con- 
tinuum physics. However, in the full electroweak theory an important ingredient 
is the Higgs field. It couples to the fermions and is responsible for particle masses 
in the broken phase. At first sight it is maybe surprising that the level crossing 
picture is still true in this case since the eigenstates of the full Dirac Hamiltonian 
are massive when one includes the coupling to the Higgs field. This is in particular 
true if we consider the ground state configuration in the broken phase : <fi = <po, 



Afj, = 0. The situation is different in the presence of a sphaleron. For a sphaleron 
configuration there is precisely one normalizable eigenmode and it has zero energy 
0,^]. As a consequence a gauge-Higgs configuration which changes continuously 
from one vacuum configuration to a neighboring one, passing through a sphaleron 
configuration on the way, will trigger an adiabatic change of the lowest positive 
eigenvalue such that it passes through zero and ends up as the "highest" negative 
eigenvalue |5|]. 

In many ways the situation around a sphaleron configuration is similar to the 
situation around the simplest non-abelian monopole. In that case there exists an 
index theorem for open three-dimensional space || which tells us that the difference 
n between the number of normalizable eigenf unctions of positive chirality and the 
normalizable eigenf unctions of negative chirality, of the Dirac Hamiltonian with 
eigenvalue zero, is related to the charge of the monopole Q by 

o = ^. (i) 

e 
n can be expressed in terms of the Higgs field: 

n=^-f (fx^ue^rd^d^, (2) 

t: /I i/Ooo 

where <p a = <p a /\4>\ for a Higgs field in the adjoint representation. Eq. (^) reflects 
that |0| is assumed to be different from zero at spatial infinity and consequently that 
the map x G S^ —>■ (f)(x) G S 2 has winding number n. In the case of the sphaleron 
we do not have a similar topological interpretation since the Higgs field is a complex 
doublet, i.e. it has four real components and we will have a map x G S 2 ^ — > 4>(x) G 
S 3 . This map is always contractable to the trivial map (^(S' 3 ) = 0). However, for 
a continuous family of static fields 0(x; r), Ai(x; r), r G [0,1], which interpolates 
between two neighboring vacua when r changes from to 1, there has to be at least 
one value r of r such that the Higgs field 0(x, r ) has a zero. In the case of the 
electroweak theory and the so-called minimal path which connects two neighboring 
vacua the sphaleron is precisely such a configuration. If we follow the minimal path 
we know that at the sphaleron there will be a normalizable eigenmode of the Dirac 
operator and it has zero eigenvalue. It is intuitively reasonable that a zero of the 
Higgs field is a necessary condition for having a zero eigenvalue of the full Dirac 
operator. However, there is to our knowledge no topological theorem which tells us 
that the full Dirac operator has eigenvalue zero for some tq. This is in contrast to 
the case of the massless Dirac operator, which only couples to the gauge field, and 
where it follows from the Atiyah-Patodi-Singer index theorem that the number of 
eigenvalues of the Dirac Hamiltonian which crosses zero is directly related to the 
change in Cher n- Simons number |7|, [|, |9|. 



Recently a number of articles have addressed the question of level crossing for 
the full Dirac Hamiltonian in the electroweak theory ]TD], |TT|, [Tj|. The proof of 
the validity of the level crossing picture has either been of numerical nature, for 
a specific chosen path in (0, At) configuration space, or it has used some special 
symmetry properties of the sphaleron solutions and is only valid for some range of 
the coupling constants. In this article we make an attempt to investigate the validity 
of the level crossing picture in the case of fluctuations of the gauge and especially 
the Higgs fields which are typical for the high temperatures of the early universe. 

There is one more reason to focus on the Higgs field. The role of this field in 
the high temperature fluctuations, and in particular its relation to the masses of the 
fermions and therefore to the level shifting, is not very clear. If the symmetry is 
restored the classical expectation value of the Higgs field is equal to zero, even if we 
have imposed complete gauge fixing. However, the field still couples to the fermions 
and the statement that (<f) T = refers to the thermal average (-) T . When we discuss 
temperature induced transitions of gauge and Higgs fields between neighboring vacua 
connected by large gauge transformations we deal with dynamical processes caused 
by thermal fluctuations, and it is strictly speaking misleading to refer to thermal 
averages. In fact we want in the discussion of the anomaly and level shifting to use an 
approximation where the gauge and Higgs fields are treated as background fields and 
even above the symmetry restoration temperature we can talk about "sphaleron-like" 
transitions. In the simulations we perform the expectation value (|^| 2 ) will typically 
be quite small compared to the classical vacuum expectation value (|y?| 2 ) at zero 
temperature, due to temperature fluctuations. This indicates that the system is in 
the symmetric phase rather than in the broken phase. It is nevertheless perfectly 
possible (for a finite volume) to talk about the fields as located in the neighborhood of 
a definite vacuum and the transition rate between such vacua. The picture becomes 
especially clear if one talks about the family of all fields which interpolate between 
the neighborhoods of two different vacua, and then focuses on those configurations 
in which the high energy thermal fluctuations have been stripped off. By such a 
procedure it becomes clear that the low frequency modes of the gauge fields can 
be viewed as interpolating between the different vacua and that the Chern-Simons 
number will change by one unit during the interpolation [Q. In this article we will 
use the spectral flow of the Dirac operator as a means to include the Higgs field in 
the analysis. 

The rest of this article is organized as follows: In sec. |2] we define the "reduced" 
fermionic part of the electroweak model we will use. In sect. || we show how to dis- 
cretize the corresponding Dirac Hamiltonian on a three-dimensional spatial lattice 
and especially how to define the concept of chirality on this lattice. Sec. f| contains 



the numerical analysis of level crossing and chirality change for a number of "hand 
made" families of gauge- Higgs configurations. This analysis is made to compensate 
for the lack of mathematically rigorous results for Dirac Hamiltonians coupled to 
chiral fermions, results which we need when interpreting the level shifting and chi- 
rality change in the numerical simulations. In sect. |5] we describe the numerical 
method used to simulate the gauge-Higgs part of the electroweak theory. Sect. || re- 
ports on the measurements of spectral flow of the Dirac operator for "sphaleron-like" 
transitions. Finally sect. [I] contains a discussion of the results obtained. 

2 The model 

Let us consider the electroweak theory where for simplicity we ignore the hyper- 
charge sector. It is of no importance for the electroweak anomaly. We also restrict 
ourselves to two families. The generalization to the twelve existing doublets is trivial. 
The fermionic part of the Lagrangian is given by: 

-$ L {h u 0XuR + h d (fXdR) - {KXur¥ ] + h d XdR^)ipL 

-$' l(K<PX!uR + KfXdR) - (Kx'uR^ + Kx'dRf^L (3) 

where the left handed doublets are 

fe = ( t: ) a -(5) (4> 



and the four right handed singlets are XuR, XdR, x'uRi Xm- With an abuse of notation 
we collect them as vectors, like the doublets, even if they of course behave like singlets 
under electroweak SU(2) transformations: 

*■=(£), M2L (5) 

The Higgs doublet tp and its conjugated field (p (which transforms as ip under elec- 
troweak transformations) are defined as usual: 

<p = (4)='^* < 7) 

There are four independent Yukawa couplings h u , hd, h' u and h' d . As mentioned 
above we have only two independent couplings in the corresponding lattice model 



and for simplicity we therefore choose h' u = h u and h' d = hd- h u ~ hd mimics the the 
quark sector (as the notation indicates), while hd = describes the lepton sector. 
It is well known that the chiral theory defined by (|||) can be given a formulation 



as a vector-like theory of Dirac fermions |I5 |. This is due to the fact that the 
representations of SU{2) are real. We can therefore introduce the charge-conjugated 
fields 

0S = £CiP' L 






where e = ir 2 acts on the isotopic indices and C is the charge-conjugation matrix, 
acting on spinor indices: ip c = Cip. In a similar way we define 



x'l = C X n= X "% • (9) 




If we use the chiral representation of the 7 matrices a Dirac fermion can be written 
as 

Id = ( { L V (10) 



and we can now introduce the two Dirac spinors: 

This allows us finally to introduce the eig/rf-component spinor \1/ by: 

.-(£). (12 ) 

and using the matrix 

M=*( h * ?), .=f « t + V (13) 



;^ 

we can write the Lagrangian ([3D as follows 



-^=«i-St T')* (14) 



such that each entry of the matrix in eq. flT4|) is a 2 x 2 matrix in isospinor indices 
and a 4 x 4 matrix in spinor indices. 

Note that flT4| ) resembles the ordinary Dirac equation but that (i) the mass is 
spatially varying, depending on if, (ii) the mass term may multiply the 1st and 
2nd (and 3rd and fourth) components differently if h u 7^ hd and, most importantly, 
(in) only the first 4 components (iPd) of the eight component spinor \I/ have gauge 
interactions. The last four (xd) are g au g e singlets. 



The notation introduced here follows [[LI]], which also introduced the generalized 
chirality matrix T 5 acting on the eight component fermion field \1/: 

f. fa o 



The matrix T^ anticommutes with the matrix in ( KM and also with the matrix 

which enters in the definition of \l/. Consequently the Lagrangian ( P| ) is invariant 
under the generalized chiral transformation |llj 



^ ^ e ^f 5 /2^ * = e iefs / 2 ^. (17) 

The symmetry expresses nothing but the fermion number conservation of the elec- 
troweak theory. In fact the current associated with the symmetry is: 

J l = *7^f 5* = ^Dlifl^D ~ XDlnlbXD, (18) 

or, expressed in terms of the original fields of the model: 

Jl = ^lI^l + ^'lIv&'l + XrI^Xr + X'rI^Xr (19) 

i.e. the fermionic current. As remarked in the introduction this current is known to 
be anomalous: 

9,4 = ^ 2 F; V F^ (20) 

and leads to the famous non-conservation of baryon number in the electroweak 
theory. The reason we introduce the generalized chirality dictated by T 5 is that 
the analogy to ordinary abelian symmetry and the level shifting picture will be 
more clear, as emphasized recently by several authors |TT], O, who proved that 



under certain assumptions a transition from one vacuum to another, related by 
a large gauge transformation, will result in a level shift. It is our goal to try to 
verify this level shifting for realistic high temperature configurations which do not 
necessarily satisfy the assumptions used in the analytical proofs of the level shifting 
and, in addition, to provide further evidence that we really observe continuum- like 
configurations during the simulations. 

3 Lattice Formulation 

First consider the time-independent continuum Dirac equation in a background tem- 
poral gauge-fixed field A,i = A^r a : 

(iaidi + Ajttj + f3m)ip(x) = Eip(x). (21) 



Here the index i runs over the three spatial directions, the wavefunction ip is a 
gauge doublet (weak isospinor) with four spinor components acted on by the 4x4 
hermitian Dirac matrices a,i and (3. These satisfy 

{a u a,} = 28a ; {a u f3} = ; f3 2 = 1. (22) 

For the purposes of this section it is again useful to consider the chiral representation 
of the 7-matrices, which means that the a,{ and /3 matrices are represented by: 

<*=(* -J;"=L, l2x2 )- < 23 > 

where 0\ are the Pauli matrices. 

On the lattice we will be using the staggered formulation for fermion fields. 

3.1 The massless case 

Let us in this section consider the massless limit m = 0. The lattice transcription 
of (|2lD reads 

i 3 

-^2m{x)[Ui{x) X {x + i)- Uj{x - i) x (x - I)} = E x (x), (24) 

1 i=\ 

where x( x ) is a single spin component isospinor field, the U{ are SU(2) valued link 
variables defined on the lattice links, which may be parameterized as exp(iaAi(x)) , 
where a is the lattice spacing, and the rji(x) are the Kawamoto-Smit phases defined as 
rji(x) = (— l) Xl ^ ^-i. It is well known that in the limit where the gauge field varies 
smoothly on the scale of the lattice spacing, the continuum limit of (23) recovers 
two copies of (plD , which we interpret as describing two independent fermion species 
15) , |i~6| . These species become coupled by terms of higher order in a |17| (i.e. by 



momentum modes from the outer half of the Brillouin zone ||18||). 

For an arbitrary background link field there are important symmetries in the 
spectrum. If we denote by (%, E) the existence of an eigenmode xi x ) with eigenvalue 
E, then 

( X ,E)^((-ir +x ^ X ,-E), (25) 

and 

( X ,E)^(r 2X *,-E), (26) 

where the Pauli matrix r 2 acts in weak isospace. Clearly applying both symmetries 
(p5|) and (^) implies that each mode is doubly degenerate. Compare with the 
symmetries of the continuum equation (|2l|) : 



(V,£)o(7W>,£), (27) 



where 75 = iaia 2 a3, and 

ty,E)<*(CT 2 r,E), (28) 

where C = ia±a^ in the particular Dirac matrix representation (|23|) . There is a 
third symmetry, which has no analogue for the lattice equation (p4|): 



(i;,E)^(pa 1 a 2 a 3 ij,-E). (29) 

Symmetry (|27|) holds only for m = 0, whereas ([28]) and (|29"D hold for any m. 

The matrix 75 defines the chirality operator in the continuum formulation. In 
order to represent chirality in the lattice system, it is necessary to project onto the 



flavor singlet sector of the multiple-species algebra implicit in (|24[) ; ie. we need to 
define an operator T^ = 75 <g> I2, where the second matrix operates on a "flavor 
space" [[17], |HJ. For the four- dimensional Euclidean Dirac operator, this projection 
has been implemented by Smit and Vink in their programme to compute topological 
charge []B5[, and by Hands and Teper in studies of chiral symmetry breaking pp[ . 
The basic algebraic steps are given in |T?|, [18| . For our three dimensional problem 
the algebra is very similar, and here we give simply the result: 

(Xn|r 5 |Xn) = ^E E (-l) X2 xi(x)U U1 Xn(x + tl+j2 + k3). (30) 

Z x ijk=±l 

Here Xn(x) is the nth eigenmode, and U\\\ denotes the average of the SU(2) link 
products over the six equivalent paths joining the two sites, thus ensuring the ex- 
pectation value is gauge invariant. Notice there is an important difference with the 
continuum operator. Since T 5 is a three-link operator, Y\ ^ 1; only in the deep 
continuum limit will T 5 and 75 coincide. At non-zero lattice spacing T 5 requires in 
principle a finite multiplicative renormalization, which may be estimated perturba- 
tively [[H|] or numerically |19| . 



Now consider what happens at a level crossing, when an eigenmode crosses zero. 
There are never any exact zero modes in a lattice simulation, so we will refer to 
modes which approach zero as "quasi- zero modes". Using the symmetries ([25]), (|26D 
and the definition (|30"D, we can see that if there is such a mode with energy e and 
chirality 7, then there will also be a mode with energy — e and chirality —7. As we 
shall see in subsequent sections, a lattice level crossing has the following signal - a 
mode initially has energy e and chirality 7 a ; as the simulation proceeds e becomes 
smaller, eventually reaching a minimum before increasing again. The chirality 7 
should reverse sign at this point, eventually reaching a value —75. Since each lattice 
mode is doubly degenerate, we interpret this process as two positive energy modes 
of chirality +7 a crossing zero to become negative energy modes of chirality +76, and, 



using symmetries (p5|), (p6|), two negative energy modes of chirality — 7 a crossing 



zero to become positive energy modes with chirality — 7&. The total change in chiral 
charge is thus AQ = — 2('~f a + 7^) In the continuum limit as T 5 approximates 75 we 
should find I7I approaching 1. Thus we find |AQ| ~ 4, which should be compared 
with the result for level crossing for the continuum Dirac Hamiltonian of \AQ\ = 2: 
the extra factor of two is a consequence of species doubling for lattice fermions. We 
shall see in the following sections how well these ideal criteria are met in actual 
lattice simulations. 

3.2 Masses and Yukawa coupling 

The situation becomes more complicated once a fermion mass is introduced. As ar- 
gued in |2~^1 , the algebra (|22"1) obeyed by a; and (3 is identical to that of the hermitian 



Euclidean matrices 7^ in a covariant formulation, which requires a four- dimensional 
lattice in order to be represented faithfully by staggered fermions. Therefore to 
incorporate mass we introduce a second "timeslice" of link variables, identical to 
the first, and couple it to the original one as if it were displaced by one unit in a 
fictitious Euclidean time direction. The lattice equation becomes 

7;J2ili(x)[Ui(x)x(x + i) - UJ (x - i)x(x - i)\ + mr] 4 (x)x(x + 4) = Ex(x), (31) 

ie. as if it were formulated on a N 3 x 2 lattice with periodic boundary conditions 
assumed in the 4th direction. The phase factor rj^x) = (— i) x i+ x 2+z 3 _ Whether 
we regard this extra degree of freedom as an extra timeslice connected to the first 
by unit gauge connections, or simply as an extra component of the wavefunction x 
is a matter of taste: the result is that we have a system which in the continuum 
limit now describes four species of four-component fermions. There are the following 
symmetries in the spectrum: 

( X ,E) & ((-i)*i+* 2 +*3+*4 X) _£ ); (32 ) 

( X ,E) «- (T 2 (-ir x *,-£). (33) 

As before, these symmetries ensure that each lattice mode is two-fold degenerate. 
However, the fact that the modes with energy —E have minus the chirality of the 
mode with energy E remains unaltered. It is also important to note that the chi- 
rality operator (|30|) links lattice sites in the same timeslice: therefore chirality for 
eigenmodes of the massive equation (|3T| ) is simply the sum of chiralities for each of 
the timeslices considered individually: 

(Xn|r 5 |Xn) = (Xln|r 5 |Xln) + (X2n|r 5 |X2n) • (34) 
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In this equation Xin(x) is an eigenvector of (|3lD projected onto the first timeslice, 
and X2n(x) the same eigenvector projected onto the second timeslice. We shall see 
that once mass is introduced, the minimum value of \E\ is m, and that the chirality 
as measured by flSTf ) tends to have opposite signs on each of the two timeslices, thus 
reducing the overall expectation. This is in accordance with our continuum intuition 
that a mass term couples fields of opposite chirality. 

Finally we discuss the incorporation of Yukawa interactions with a Higgs field, so 
that the standard model can be discussed. Once again we begin with the continuum 
formulation: if we use the explicit representation of the Dirac matrices fl2"3|) we may 
according to the Lagrangian ([14]) write the equations of motion (suppressing the 
primed fields, which are not essential to the formulation) as 



E r, L , (35) 



M f -ioidi) \xr) " \Xr 

where each entry is a 2 x 2 matrix in spinor space and a 2 x 2 matrix in isospace. 
The resulting system is very similar to the Dirac equation for /ottr-component Dirac 
spinors (recall the comments following (|T4D). Now, to model a mass term in the 
original system, i.e. to incorporate a prescription for the matrix (3, recall that we 
doubled the number of spinor components by introducing a second timeslice with 
identical gauge interactions, and had the mass terms couple the two slices. Here 
we will do the same, except that on the second timeslice, corresponding to the \ 
components, there will be no gauge interactions. Gauge invariance will be enforced 
as in the continuum by coupling the Higgs fields to the ip fields on the original 
timeslice. The system of lattice equations then reads 

I Ei rn{x)[Ui(x)il>(x + i)- UJ(x - i)ip(x - i)} + r] 4 (x)M(x) X (x) = Eif>(x), (36) 
EiVi( x )[x(x + i) -x(x-i)] +r) A (x)M^(x)ip(x) = Ex(x), (37) 



2 



where if) and x are now single spin- component gauge-doublet fields. The following 
symmetries exist: 

U>, X ,E)&(ri4,iJ>,-n4X,-E), (38) 

and, if h u = h d , 

(iJ;,X,E)&(T2il>*,-T2X*,-E)' (39) 

Once again, these symmetries ensure that each lattice mode is doubly degenerate if 
h u = hd', if not, there will be an additional fine structure in the spectrum. 

Because of the symmetry of the spectrum about E — 0, the system described 
by the lattice Hamiltonian (J36l)- (|37|) is not truly chiral: if it were there would be 
spectral asymmetry, and true spectral flow would result when a mode crossed zero. 

11 



Of course, the reason the model is not chiral is due to the species replication inherent 
to the lattice approach (we leave aside the issue of the requirement for complex 
representations, since in any case it is perfectly possible to introduce a U(l) y gauge 
interaction in the lattice system). For every positive energy level which crosses zero 
to become negative, there will be a corresponding negative energy level crossing zero 
to become positive. The symmetries (|38|) , (|39|) ensure that at each zero crossing four 
SU(2) doublets are involved: two cross one way and two the other. Thus there is no 
true baryogenesis on the lattice, since for each unit of baryon number created there 
is another one destroyed. Of course, in the continuum limit we expect the species 
to become independent (in the quenched approximation), and thus we should be 
able to focus simply on the crossing itself as the event of interest. Note also that 
chirality is still defined by (0), and that we have the freedom to either define the 
total chirality of the mode as {ip\r 5 \ip) + (x\T 5 \x) , or the more sophisticated version 
(f/>|r 5 |'?/>) — (x|r 5 |x). As we saw in sect. ^ it is the last version which is directly 
related to the anomalous fermion number of the electroweak theory and we will use 
this definition in the following: 

*|f 5 |*) = (VIW) - ( x |r 5 |x) , * = (^) (40) 

4 Calibration 

As mentioned above there seem not to be any analytical results which tell us in 
detail about the level crossing for chiral fermions coupled to the gauge-Higgs system 
encountered in the standard model. In order to know what to expect we have 
investigated this question numerically on the lattice by taking families of smooth 
configurations^ which interpolate between different classical vacua and for these 
configurations we have followed the flow of eigenvalues and the change of chiralities. 
We have investigated the following scenarios: 

(1) A family of gauge-Higgs field configurations which interpolate between two 
vacua and where the zero of the Higgs field appears at the same time as the 
Chern-Simons number N cs (t) of the gauge field is equal to 1/2. 

(2) A family of field configurations where the gauge- and Higgs fields still inter- 
polate between two vacua, but where the zero of the Higgs field is displaced 
in time relative to the time t where N cs = 1/2. 



4 By "smooth" configurations we mean lattice configurations which are approximations to con- 
tinuum configurations. 
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(3) A family of field configurations where only the gauge field interpolates between 
two neighboring gauge vacua, while the Higgs field remains constant. 

(4) A family of field configurations where only the Higgs field interpolates between 
two neighboring Higgs vacua, while the gauge field stays constant (for instance 
zero). 

The field configurations are constructed as follows: Let us first define a gauge trans- 
formation V(x) with winding number one: 

V(x) = (-l)exp (|^^| m ax{|^(x)|,j = 1,2,3}) , (41) 

where x denotes a lattice point, L the linear size of the lattice, Vi = Xi — \L and oi 
are the Pauli matrices. The trivial vacuum is given by 

Ui°\x) = i, (O) (*)=(J), (42) 

while a vacuum configuration with winding number one is given by 

U\ l \x) = V(x)V- 1 (x + e i ) = exp(ia i u i (x)) (43) 

^V ( x ) = V(x)(/)^. (44) 

where e^ is a unit lattice vector in direction % and Ui(x) will be close to zero for large 
lattices. 

A continuous family of configurations which interpolate between the two vacua 
is simply given by 

U i (x,t)=exp(if(t)cr i Ui(x)), (45) 

<P(x,t) = (l-g(t))^ + f(t) ( f ) (1 \x), (46) 

where /(£) and g(t) are monotonic, continuous functions which interpolate between 
and 1 when t changes from to the final time which we denote T. 

Note that there exists precisely one space-time point (x, t) where <j)(x, t) = and 
the time t is determined by 1 — g(t) = f(t). We denote this time th- If g{t) = f(t) 
it will coincide with the time t\/2 where N cs (t) = 1/2, but if g ^ f it will in general 
be different from t\ii- 

We will choose f(t) = t/T in the actual numerical calculations. This means that 
ti/2 = T/2. In case (1) and (4) mentioned above we choose g(t) = f(t). In case (2) 
we choose g(t) = min{3t/T, 1}, which means that the zero of the Higgs field will be 
located at t — T/4, (t h = T/4) i.e. displaced by T/4 relative to ii/ 2 . 

13 



Let us now describe the spectral flow and the change in chirality for the four 
families of field configurations mentioned above. In fig. [I] we show the Chern- 
Simons number as a function of the time t for the gauge fields given by (|45|). In 
addition we show in fig. ^ the change in the lowest eigenvalue and the change in 
chirality for various values of the Yukawa coupling. Everything behaves as expected 
from the continuum considerations. 

In case (2) where the zero of the Higgs field is displaced from ti/2 we see that 
the time to at which the lowest eigenvalue crosses zero changes continuously with 
increasing Yukawa coupling from tin, corresponding to the crossing at N cs = 1/2 
for zero Yukawa coupling, towards t^, which is the time where the Higgs field has 
a zero. This is illustrated in fig. |^. The flip in chirality always occurs at t , as 
expected. 

In case (3) where the Higgs field is constant we see that the level crossing is still 
present for moderate values of the Yukawa coupling. This is in sharp contrast to 
the situation where we insert an explicit mass term in the Dirac equation. In that 
case we know from the continuum equation that there is a mass gap, separating the 
negative and the positive eigenvalues, and this is also seen explicitly in the lattice 
implementation (see sec. H). When the Yukawa coupling is increased t is gradually 
displaced towards T, and for a finite value of the Yukawa coupling t reaches T and 
there is no longer any level crossing and no chirality flip. This scenario is shown in 

fig- I 

The behavior of the lowest eigenvalue in the situation where the gauge field is zero 

(i.e. Ui(x) is given by fl4"2])) and Higgs field is given by fl4"6| ) with g(t) = f(t) = t/T 
is shown in fig. |^ as a function of the Yukawa coupling. The situation is in a 
certain sense dual to the one in case (3). For small Yukawa coupling there is no 
level crossing, but eventually a zero and a corresponding chirality flip takes place at 
T and with increasing Yukawa coupling t moves towards t h . 

The results in this section are relevant when we want to interpret the data 
obtained in real lattice simulations. In particular we see that increasing the Yukawa 
coupling can displace the time to of zero level crossing. It is somewhat surprising 
that there is no need for a zero in the Higgs field in order to observe level crossing, 
provided the Yukawa coupling is sufficiently small. However, our three dimensional 
lattice has periodic boundary conditions and since it is finite, it imitates a compact 
space of finite volume. This means that the "continuum" configuration imitated by 
(f45|) for the gauge field and fl4~2"D for the Higgs field will have finite energy, and we 
expect spectral flow of the continuum Dirac operator to be continuous with respect 



14 



to the coupling constants in the lagrangian. Since we have level crossing for zero 
Yukawa coupling by the Atiah-Singer index theorem and since the chiral symmetry is 
not broken for finite Yukawa coupling, we must have level crossing for small Yukawa 
couplings as well, even if the Higgs field is nowhere zero. This is indeed what we 
observe. For sufficiently large Yukawa coupling the level crossing can disappear 
without violating any continuity requirements by the mechanism shown in fig. |j. 
It is more difficult to decide whether these continuity arguments remain valid if the 
three dimensional space is non-compact. The very simple configurations we have 
used will have infinite energy in the non-compact case and might not belong to 
a reasonable class of field configurations where one can define a spectral flow of 
the Dirac operator H (Ai(t) , <f>(t)) . Similar remarks are valid in case (4). On the 
other hand the results obtained in case (2) are probably true also for a non-compact 
continuum space, even if if the actual configuration used here has infinite energy in 
the continuum limit for t 7^ 0,T. The reason is that it is easy to find finite energy 
continuum field configurations with the characteristics used in case (2). An analytic 
proof of the behavior observed in case (2) has recently been given if space is one 
dimensional ||26|| . It would be interesting to extend it to three dimensions. 

At high temperature we might certainly encounter situations where the phase 
of the Higgs and the gauge field are placed in different gauge sectors as in case 
(3) or (4), since the finite temperature provides the needed infinite energy. Conse- 
quently cases (3) and (4) are relevant for our interpretation of data not only because 
the spatial lattice used in the simulations is finite, but also because the situations 
might be generic at the high temperatures prevailing in the early universe above the 
electroweak phase transition. 

5 The numerical method 

5.1 The approximation 

The following approximations have been used in this work: we consider only one 
weak isospinor fermion and the two associated singlet fermions and we ignore the 
hypercharge sector. These approximations are presumably not important for the 
questions we want to address. In addition we completely ignore the feedback of 
fermions on the gauge and Higgs fields and we always treat the gauge and Higgs 
fields as background fields when we ask dynamical questions about the fermions. 
This approximation is difficult to control, but it might still not be very important 
for the dynamics of the gauge and Higgs fields. The basic approximation, when it 
comes to the dynamics of the gauge and Higgs fields, is the use of classical physics 



15 



as the tool for simulating the high temperature fluctuations in the electroweak the- 
ory. One would expect this approximation to give a reliable representation of the 
magnetic sector of the theory since the classical partition function on the lattice^ 
is identical to the naive infinite temperature limit of the full quantum partition 
function. If we discuss the electric properties of the theory it is by now well known 
that we need to supplement the naive dimensional reduction with mass and charge 
renormalizations and in addition with Higgs fields in the adjoint representation (the 
remnants of the Aq component) in order to get a correct representation of the full 
quantum theory Jl4| . The Gauss constraint will be satisfied by the use of the classical 
equations, and this means that the kinetic part of the adjoint Higgs field coupled to 
the gauge field will be included in the classical approximation (The Gauss constraint 
is equivalent to a term (D i Ao) a (DiAQ) a in the Lagrangian). However, there will be 
additional terms like ^4q'^o an d ^o i n the effective three-dimensional high temper- 
ature expansion of the action, where especially the first two terms are important 
for the Debye screening mass. We will ignore these subtleties here[] and view the 
present computer simulations as generating typical high temperature fluctuations, 
the temperature being judged to correspond to fluctuations in the unbroken phase 
since ((p 2 ) is significantly smaller than its tree value for the choice of lattice coupling 
constants we will use (see later). 

Classical thermal fluctuations can be generated either by using a canonical en- 
semble and applying some standard Monte Carlo updating like heat bath or the 
Metropolis algorithm to the field configuration, or (more fundamentally) to use a 
micro-canonical ensemble and just let it develop according to the classical equations 
of motion. In the first case the temperature is kept fixed, while the energy per de- 
gree of freedom will have fluctuations proportional to 1/yN, N being the number 
of degrees of freedom. In the second case the energy will be constant, while the 
temperature of the system will fluctuate as 1/y/N. From a numerical point of view 
the first method is preferable, since the second is known to suffer from problems with 
ergodicity for too small systems. We have nevertheless chosen to use the classical 
equations of motion since they reflect very directly the real time evolution of the field 
configurations. However, in order to use the classical equations of motion we have to 
prepare an initial distribution of fields and their conjugate momenta corresponding 

5 Let us remind the reader that there exists no classical partition function of a continuum field 
theory due to the Raylcigh-Jeans instability. If one restricts the number of degrees of freedom per 
unit volume by considering a lattice version of the field theory the concept of a classical partition 
function is perfectly well defined. 

6 The role of these additional finite temperature corrections is not clear for the kind of questions 
we will deal with. Leaving aside the question of a gauge invariant meaning of the Debye screening 
mass, it is a static quantity, calculated as a temperature average in a canonical ensemble, while we 
consider dynamical processes such as sphaleron transitions which occur in real time. 
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to a given temperature, and here we have to use the heat bath or the Metropolis 
algorithm. For systems with gauge invariance this means additional complications 
since the Gauss constraint has to be satisfied for the classical configurations and 
this will not be the case for a Metropolis updating. This technical obstacle can be 
circumvented by introducing the Gauss constraint as a Lagrange multiplier (we refer 
to for details). After the initial thermalization by means of Metropolis, we let the 
system develop according to the classical equations, discretized in a way which auto- 
matically conserves the Gauss constraint (see next subsection). For sufficiently large 
systems the fundamental hypothesis of ergodicity of the micro-canonical ensemble 
will ensure that the phase space surface of constant energy is covered uniformly if 
the system evolves according to the classical equations of motion. 

5.2 The gauge-Higgs system 

As mentioned above the sequence of gauge and Higgs field configurations we are using 
is generated by the classical equations of motion. By that we mean the following: we 
choose the standard Wilson action for Euclidean lattice gauge theory combined with 
the standard Higgs action on a lattice and make the obvious change of sign of the 
spatial part relative to the temporal part in order to get a Minkowskian signature. 
After that we make the lattice spacing in the temporal direction much smaller than 
in the spatial directions by scaling the lattice spacing "a" to "aAt" . In this way the 
action will be split into a kinetic and a potential part: 

S = J2( E kin - E pot ) (47) 

t 

where the summation is over all time slices. From S one can derive the classical 
equations of motion which connect timeslice t and t + At. These equations of motion 
will automatically satisfy the Gauss constraint. 
If we define an energy functional H by 

H = E kin + E pot (48) 

it will not be strictly conserved by the equations of motion since time is still discrete. 
However, in the limit At — > 0, H reduces to the correct Hamiltonian of the system 
and the energy will be conserved by the equations of motion. For At < 0.05 we 
have, for the range of coupling constants used in the simulations, found very good 
conservation of the energy. 

For a detailed discussion of the equation of motion derived from (|47|) we refer 
to ; let us only here give explicitly the form of J2t Ekin and J2t E po t used in the 
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simulations: 



£3 



pot 



J^-y £ Tr (*t ^ _ $t ^ ^ +0$x+e 

n(s) z z i,t,i 



(49) 






(50) 



where U XjX+fJ , is the lattice gauge field on the link connecting x and x + fi, U n denotes 
the plaquette action for the gauge field and □ (£) stands for a plaquette in the — % 
plane while C(s) refers to a plaquette in a % — j plane. The SU(2) Higgs doublet tp 
is represented as a matrix 

H4&)- (51) 

The summation over spatial points is finite and limited by the lattice 3-volume 
(periodic boundary conditions) while it is infinite in the time direction. 
The Metropolis updating used to thermalize the system uses as weight 



exp(-i^), H^ = H + iG 2 



(52) 



where H is the lattice hamiltonian given by ([§), while G denotes the lattice version 
of Gauss constraint (see @ for details). By choosing £ sufficiently large one can (at 
the expense of slowing down the thermalization time) control the violation of Gauss 
constraint to be arbitrarily small in the Metropolis updating. 
If we write the classical continuum Hamiltonian as 



H 



(TX 



l -E?Et + \F^ + M 2 + | DM 2 + A(M 2 - v, 



2\2 



the tree value connection between the lattice parameters in (p9|). 
tinuum coupling constants in (|53|) is as follows[] 



M 



2(1 - 2[3 R - 3[3 H ) 
f3 H a? 



X, 






i: 



4 



(53) 
and the con- 

(54) 



where a denotes the lattice spacing and 



2 2 

t> a 



Ph 



Ph_ 



(2fo + 3/3 H -l). 



(55) 



7 We keep this unnecessary complicated notation for historical reasons 
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The naive continuum limit a — > imposes a fine tuning on {2(3 r + 3(3h)- If we 
ignore thermal fluctuations and determine the jy-mass My/ and the Higgs mass Mh 
from the coupling constants in the Hamiltonian we have: 



,, 2 _ 16(3 R 2 2 2 _ @H2 



M 2 H a 2 = ^Pv\ M^a 1 = ^V. (56) 

Ph Pg 

If for convenience, to restrict the coupling constant space, we impose the constraint 
M w = Mjj, we get the relation 

'« = isfc < 57) 

From H we know that a reasonable choice is (3q > 10 and (3h in the neighborhood 
of 1/3. This means that (3r is quite small. 

5.3 The procedure 

Equipped with an initial configuration at a given temperature and the classical 
equations of motion we can now follow the evolution of the Chern-Simons number 
by measuring the integral 

Q{t) = N cs {t) - N cs {0) = JL jf * dt J d 3 x F^F**. (58) 

This quantity has a simple implementation on the lattice and will directly give us 
the change in Chern-Simons number, which is a gauge invariant quantity. This 
procedure has been applied succesfully before (see p| for details). The typical mea- 
surement of Q{t) will result in a curve which stays for some time around integer 
values of Q{t), while the transitions between these values take place rather rapidly. 
Superimposed on these movements of Q{t) will be short wavelength thermal fluctu- 
ations, which occasionally can mask somewhat the picture of Q{t) jumping between 
integer values. If desired, it is possible to strip off these thermal fluctuations in 
the following way. By solving the classical equations of motions we get a time se- 
quence of configurations. For each of these configurations we can apply the simplest 

relaxation equation: 

dip _5H_ dA _5H_ 

~dt ~ ~ V ~dt ~ ~Ja { ' 

where H is the Hamiltonian (H8). This technique is well known from the lattice study 



of monopoles and instantons |20| , p2| , |23| , p4] , p5| . If one "cools" the configurations 
too much one forces them to a vacuum configuration, which is not what we want. 
However, only a few cooling steps will reduce the energy of the gauge configurations 
by a factor 80 without breaking the link to the original configuration since the 
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part of the field configurations which survives this relaxation is precisely the long 
wavelength fluctuations which we imagine are responsible for the large scale changes 
of Q(t). Indeed, the picture of the transitions between different integer values of 
Q(t) is considerably sharpened by such a cooling. This effect is illustrated in fig. ||. 

In order to check for level crossing during the time-development of the gauge 
field configurations we calculate the eigenvalues of the Dirac Hamiltonian for each 
time step, using the staggered fermion formalism discussed in detail in sections @ 
and |3[ Since the Dirac operator Hr)(A(t), (p(t)) will depend implicitly on the time t 
due to the time dependence of the gauge and Higgs field we will get a spectral flow 
of the eigenvalues and we can compare this with the expectation that a crossing of 
zero should be related to a change in Chern-Simons number and that diving and 
rising eigenvalues have different generalized chiralities. 

The spectrum of the Dirac operator was found numerically by first using the 
Lanczos algorithm to tridiagonalize the hermitian Hamiltonian matrix, and then 



using Sturm bisection to extract the eigenvalues. As described in [20], on large sys- 
tems this procedure finds eigenvalues iteratively, with convergence tending to occur 
first for extremal values of \E\. We found that for the systems we investigated the 
smallest eigenvalue converged after roughly 300 iterations of the Lanczos algorithm, 
and the lowest ten eigenvalues after roughly 700 iterations. Once the eigenvalues had 
converged the eigenmodes of the tridiagonal matrix are easily extracted by inverse 
iteration: the Lanczos algorithm is then rerun with the same initial vector in order 
to perform a unitary transformation on the eigenvectors of the tridiagonal matrix 
to convert them to eigenvectors of the original Hamiltonian. We checked that the 
resulting residual vectors had norms in the range 10~ 6 - 10~ 10 , which is adequate 
for the subsequent measurements of chirality. 

6 The measurements 

6.1 Results for the massless Dirac operator 

In [0 the spectral flow of the massless Dirac operator^ was investigated along the 
lines outlined above. It was shown that there indeed was a crossing of zero associated 
with a sphaleron-like transition where the Chern-Simons number changed by one 
unit. Some questions were left unanswered, however. 

The first question was connected with the chirality. As explained in sect. |3] the 
spectrum is symmetric with respect to energy E = 0. If an eigenvalue of chirality +1 

8 By "massless" we mean that the mass is explicitly chosen zero and that there is no coupling 
to the Higgs field. 
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crosses zero from above, a corresponding eigenvalue of chirality -1 should cross zero 
from below. This is indeed what we observe. In fig. |6] we have shown the change 
in Chern-Simons number which characterizes a typical sphaleron-like transition. In 
fig. [7| we show the corresponding diving of the lowest eigenvalue and the associated 
change in chirality. 

The second question was related to the observation that a number of configura- 
tion histories had a Chern-Simons number which changed by 1/2 and then returned 
to its original value. These changes survived relaxation and they should consequently 
have some topological content. This hypothesis was supported by observing that 
their appearance was accompanied by a diving of eigenvalues to zero. It was conjec- 
tured in |3] that one could view such configurations as gauge-Higgs configurations 
which climbed to the sphaleron barrier where the Chern-Simons number is 1/2, but 
then rolled back down to the original vacuum rather than to the neighboring one. 
By measuring the chirality we can substantiate these conjectures. This is illustrated 
in fig. [8], where we show such a situation. The eigenvalues vanish twice, the two 
zeroes being quite close, while the chirality at the first zero changes from +1 to -1 
and back again to +1 at the second zero. The obvious interpretation is that the field 
configurations reach a "sphaleron" -like configuration (where the energy is zero) and 
that it also overshoots it a little causing the negative chirality state to cross zero, 
but afterwards it returns to its original vacuum in agreement with a total change of 
Chern-Simons number of zero, and once again the chirality states cross zero. 

We conclude that the gauge fields seem to qualify as continuum configurations. 
The change in Chern-Simons number is reflected closely in the spectral flow of 
the eigenvalues of the massless Dirac operator, which is implemented by means of 
staggered fermions. 

6.2 Massive fermions 

Before we turn to the "complete" electroweak theory, where the fermions acquire 
their masses via the Higgs field, it is illuminating to address the coupling of the 
gauge field to the massive fermions as discussed in sec. |3]. In this case the continuum 
fermion is a four-component Dirac spinor and has isospin 1/2. Recall that the lattice 
implementation uses two fields with support on alternating time slices to represent 
the lattice version of the Dirac matrix /3, and chirality for a given eigenmode is 
defined by (HE 

In this case the eigenvalues of the Dirac operator do not cross zero if solved in 



9 It makes no sense to use the generalized chirality T$ since the spinors in this case do not split 
into different isospin representations 
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a sequence of gauge field configurations where the Chern-Simons number changes 
by one. In fig. |S]a we show the lowest eigenvalue for a number of different masses 
and we see that the lowest eigenvalue is approximately equal to the mass m, in 
agreement with expectations from the continuum formalism. We have studied the 
chirality of the lowest mode as a function of the mass. Some results are shown in fig. 
||b. If the mass is zero the component X2n(x) (see (|3^)) on the second timeslice does 
not contribute to the chirality as defined by ([34]). However even for small masses 
both components will start to contribute, but the sum for small masses is quite 
similar to that obtained in the massless case. However, as the mass increases (T^) 
decreases towards zero, due to a cancellation between the contributions from the 
two timeslices. For the lowest eigenmode it seems still possible to observe a change 
in sign of (rs) even for large masses. If we look at higher eigenmodes the chirality 
is closer to one and is less affected by the mass term. All this is in agreement 
with expectations from the continuum theory and verifies that the "two-timeslice" 
formalism works satisfactorily. We are now ready to apply the formalism to the 
electroweak theory. 

6.3 Spectral flow of the full Dirac operator 

In the following we will assume, unless explicitly mentioned, that h u = h^ and we 
will denote the common coupling constant by h y . We have investigated the level 
crossing picture and the change in the generalized chirality T^ for a number of se- 
quences of configurations, each sequence corresponding to a change in Chern-Simons 
number of order one. For each of these sphaleron-like transitions the eigenvalues and 
eigenmodes of the Dirac equations (|36"D , fl3"T|) were found for different values of the 
Yukawa coupling h y ranging from 0.05 to 1. In addition these measurements were 
performed for three different "versions" of the same sphaleron-like transition, cor- 
responding to different degrees of relaxation of the individual configurations, more 
precisely 0, 4 and 8 relaxation steps in the discretized version of eq. (^9|), which had 
the effect of stripping off as much as 96 % of the short wavelength thermal energy 
in the case of 8 relaxation steps. 

The picture is very consistent for the various transitions and we will present 
details from two of them. In fig. |i~0| a and |i~T| a we have shown the lowest eigenvalue 
for the two sequences of configurations (parametrized by the time t) near sphaleron 
transitions for Yukawa couplings h y = 0.1, 0.25 and 0.5 (fig. |10|a) and h y = 0.1 and 
0.5 (fig. |Tlla). In fig. p~0J 3 and lib we have shown the measured generalized chirality 



(r 5 ) for the same sequences of configurations and the same values of the Yukawa 
coupling. 8 relaxation steps have been imposed on the configurations. The next to 
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lowest eigenvalues are only weakly affected and the same holds for the corresponding 
(r 5 ) values and we have chosen not to show these. 

The effect of the relaxation is to increase the value of (f^) slightly, and to affect 
only rather weakly the lowest eigenvalues in the case where h y = 0.1. If h y = 0.5 
the effect of cooling is somewhat stronger. This is show in fig. [T3 for the transition 



shown in fig. 11. In fig. |T2| a and p"2] b h y = 0.1 and three curves are shown of the 
lowest eigenvalue and (I5), respectively, corresponding to 0, 4 and 8 relaxation steps. 
The lowering of (V5) in the case of no relaxation is in accordance with the discussion 



in sec. and is presumably caused by fluctuations of the Higgs field. In fig. 12c and 



12| d we show the same set of curves, but generated for h y = 0.5. It is clear from the 
figures that it will be difficult to conclude anything about level crossing and chirality 
flip of h y > 0.5 without relaxation of the individual configurations. When h y > 1 it 
is impossible even with extensive relaxation. For such large values of h y there seems 
to have been introduced a lot of spurious transitions due to the fluctuations of the 
radial part of the Higgs field. 

It is of interest to understand the relationship between the Higgs field and the 
transition between positive and negative chirality. As we increase the Yukawa cou- 
pling we indeed observed some change in the precise location of the "zero" energy 
solutionis as is clear from fig. |l0| b and |i~l"]b. The change is small compared to 
the total width of the sphaleron-like transition (the region where the Chern-Simons 
number change by approximately one unit), but it nevertheless shows that there 
cannot be a strict relation between the zeros of the Higgs fields and the occurrence 
of a zero mode. This is in complete agreement with the analysis performed in section 
f|, where we used artificial "smooth" configurations. On the other hand we have to 
say that our setup is not the best for an investigation of such a relationship. The 
reason is that we are in the symmetric phase of the electroweak theory. The ex- 
pectation value of the Higgs field is quite small and it fluctuates significantly. This 
means that 4>(x, t) might get close to zero in a number of places and the relaxation 
we use is insufficient to change the expectation value of the Higgs field very much. 
The phases of the Higgs and the gauge field seem to align much faster. This is due 
to the small value of the Higgs coupling which is dictated by the formal continuum 
limit. It would be interesting to perform the same measurements on configurations 
where the temperature is below the transition point, such that we are still in the 
broken phase. Work in this direction is in progress. 

Finally we have investigated the situation with different Yukawa couplings for 



10 As remarked previously a "zero" energy mode should be understood as a quasi-zero mode. 
The discrete approximations used in the numerical calculations imply that we have to be content 
with eigenvalues which get (very) close to zero. 
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the upper and lower components (h u = 0.1 and hd = 0.25). In this case the Dirac 
Hamiltonian ceases to be charge conjugation symmetric, and as discussed in sec. || 
we expect states which were previously degenerate to split into two states of different 
energy. To test this we calculated the inner product < ^j | r/ 4 T 2 \I/* >, where i and j 
denote the two eigenmodes of nearly equal energy. According to the symmetries (|38D , 
(P^), the previously degenerate eigenmodes would have been related as *£fj = r/4T 2 \I/* 
in the limit h u = h^. As expected, we find the inner product to be close to unity for 
all such pairs. 

7 Discussion 

The purpose of this work has been to establish the possibility to investigate in detail 
by numerical methods a number of questions related to sphaleron-like transitions 
in the early universe. By simulations we have generated conditions which should 
have some resemblance with the conditions near the electroweak phase transition: 
a small expectation value of \ip\ 2 (compared to the tree-value determined from the 
Lagrange function), but the temperature still not so high that it completely masks 
the transitions between neighboring vacua. We are here in the lucky situation that 
we can use our finite lattice volume in a constructive way. Had we been working 
at infinite or very large volume V it is clear that the temperature fluctuations in 
globally defined quantities like the Chern-Simons number would be very large since 
they essentially grow with the square root of the volume. By fine tuning the volume 
and the coupling constants we can work at high temperature and corresponding 
small expectation values of \<p\ 2 , but still have an acceptably small change in the 
Chern-Simons number . The risk is that finite volume artifacts will be present, and 
certainly we would have to check this aspect carefully if we attempted to measure 
actual transition rates. Then one would have to make sure than the transition rate 
grew correctly with volume. In this study we have concentrated on the qualitative 
aspects of the individual transitions: is there any trace of topology left in the high 
temperature phase above the electroweak transition and can it be associated with 
level crossing in the usual way? 

We have seen that the gauge fields indeed behave like continuum configurations 
with respect to change of Chern-Simons number and level crossing and chirality 
of the massless fermions. As mentioned in the introduction it is not a priori clear 
how this picture survives the coupling to the Higgs sector of the electroweak theory. 
However, for not too large values of the Yukawa coupling h y ( h y < 0.5 ) we have 
observed the diving of a single eigenvalue to zero and an associated flip in the 
generalized chirality T 5 when the Chern-Simons number changes by one unit. The 
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natural interpretation of this is that energy eigenvalues of eigenmodes of opposite 

generalized chirality cross at zero energy. The crossing is not linked in any precise 

way to the point where ANcs = 1/2, but as discussed in sec. f|one would not expect 

this for general gauge-Higgs configurations^. The observed picture of level crossing 

and especially chirality flip is enhanced by stripping off the high energy thermal 

fluctuations of the gauge and the Higgs fields. This is an indication that the low 

frequency part of the Higgs field indeed has some topological content, even if it is 

clear that the Higgs configurations from the outset are rather "wild" . It also shows 

that we presumably have a "renormalization" of the chirality due to fluctuations, as 

mentioned in sec. |3|. The break down of any sensible eigenvalue flow and chirality 

flip for h y > 1 is probably due to the fluctuations of the radial part of the Higgs 

field. In sec. f| where the radial part of the Higgs field was smooth we could choose 

a considerably larger h y before we encounted any problems with eigenvalues and 

chirality. When we use the relaxation algorithm on the configuration the situation 

does not improve much. On the other hand we know from direct measurements 

that the phase of the Higgs field and the gauge field are almost totally correlated 

after 8 relaxation steps and that we get very sensible results for smaller values of 

h y . Since our relaxation algorithm works rather slowly on the radial part of the 

Higgs field, it is natural to expect the fluctuations of the radial part to develop into 

lattice artifacts for some threshold value of h y . This threshold value should depend 

somewhat on the relaxation and we have observed this dependence. 

Finally we wanted to use the measured flow of eigenvalues and the chirality flip 

to extract information about the topology of the Higgs field. The results obtained 

in sec. [| made this procedure more ambigious than desired. Before relaxation we 

can by direct measurement see that the phase of the Higgs field and the gauge 

field are rather decorrelated. After relaxation they are highly correlated. The basic 

picture of levelcrossing and chirality flip is not changed by relaxation for h y < 0.5 

although it becomes much clearer. Unfortunately we have seen in sec. ^ that we can 

have level crossing and chirality flip even if the phase of the gauge and the Higgs 

field develop into different vacua. This means the the cooling itself could induce a 

change of the relative phase of the gauge and the Higgs field (from being in different 

vacuum sectors to being in the same vacuum sector) without a drastic change in the 

levelcrossing and chirality flip. To determine in more detail the topology associated 

with the Higgs field would therefore require additional independent measurements. 

11 Even for Yukawa coupling zero we are not aware of any analytic considerations which tell us 
that the crossing has to take place at the point where ANcs — 1/2- Furthermore our configurations 
do not interpolate between two vacuum configurations due to the finite temperature. The best we 
can say is that we start out in some neighborhood of a classical vacuum and end up near another 
classical vacuum. 
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Figure Captions 



Figure 1: The Chern-Simons number as a function of the time t for the gauge 
configurations fl45|) with f(t) = t/T. 



Figure 2: (a): The lowest eigenvalue as a function of time t for the gauge configu- 
rations ([45]) (case (1)) h y = 0, 0.5, 1.0 and 2.0. (b): The chirality T 5 for the same 



values of h y = as in (a) 



Figure 3: (a): The lowest eigenvalue in case (2) for h y = 0, 0.5, 1.0 and 2.0 as a 
function of the time t. (b): T 5 for the same values of h y . 
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Figure 4: (a): The lowest eigenvalue in case (3) for h y = 0, 0.25, 0.5, 0.75, 1.0 and 
2.0 as a function of the time t. (b): T 5 for the same values of h y . 



Figure 5: (a): The lowest eigenvalue in case (4) for h y = 0.1, 0.25, 0.50, 0.75, 1.0 
and 2.0 as a function of time t. (b): T^ for the same values of h y . 



Figure 6: A typical "sphaleron" transition before relaxation and after 6 relaxation 
steps. 



Figure 7: (a): The "crossing" of the lowest eigenvalue and (b): the change in chirality 
for the transition shown in fig. 1. 



Figure 8: (a): The lowest eigenvalue for a sequence of configurations which reach the 
"sphaleron peak" corresponding to AN C s = 1/2, and the return to configurations 
with ANqs = 0. (b): The change in chirality for the same configurations. 



Figure 9: (a): The lowest eigenvalue Ao of the massive Dirac equation ([H|) for 
different choices of the mass close to a sphaleron-like transition. On the figure we 
shown the difference Ao — m which is always positive, (b): The chirality (|34"D for the 
same configurations and masses as in (a). 



Figure 10: (a): The lowest eigenvalue near a sphaleron-like transition for h y = 0.1, 
0.25 and 0.5. (b): (r 5 ) for the same configurations and same values of h y as in (a). 



Figure 11: (a): The lowest eigenvalue for near a sphaleron-like transition for h y = 0.1 
and 0.5. The time to where the eigenvalue crosses zero has a much weaker dependence 
on h y than the transition considered in fig. [10]. (b): (T^) for the same configurations 
and values of h y as in (a). 



Figure 12: (a)-(b): The lowest eigenvalue and (f 5 ) for the sphaleron-like transition 
used in fig. [11] for h y = 0.1, if the configurations has undergone 0, 4 and 8 relaxations. 
(c)-(d): The same figures as in (a)-(b), only with h y = 0.5. 
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Fig. 3 
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